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SUMMARY 


The main objective of this program was to develop and flight test a high 
profile 4-D descent. A profile descent is an aircraft’s descent from cruise 
altitude to the terminal area at or near flight idle, as shown in figure 1. 
The high profile descent is the upper part of the descent, between cruise and 
a metering fix, the remainder of the descent from the metering fix to the 
runway is called the runway profile descent. A 4-D high profile descent is 
one in which an aircraft begins its descent at a specific point such that it 
arrives above a point on the earth at a desired altitude - and at a specific 
time . 


Software in the existing production L-1011 Flight Management System, 
which has a 3-D capability, was modified to perform 4-D descents in flight. 
Eleven descents were performed — eight on flights off the coast of California 
between June 25 and July 30, 1979, and three on a demonstration flight to the 
Dallas/Ft. Worth Regional Airport on August 1, 1979. This last flight was 
undertaken to demonstrate the 4-D descent capability in the arrival airspace 
of a major hub airport. It was desired to exercise the system in conjunction 



Figure 1. - Descent profile. 




with the computerized time-based metering system that is used in the Ft. Worth 
Air Route Traffic Control Center to control flow to the Dallas/Ft. Worth 
Regional Airport. This system was mainly developed at the Ft. Worth Center 
and is now also used at the Denver Center. It is planned to make it available 
to all FAA centers in 1980. 

Early in the program an agreement was reached between NASA and Lockheed 
that the end-of-descent arrival time accuracy goal would be ±30 seconds. Other 
goals were: airspeed 250 ±5 kts, position ±1 n.mi., and altitude ±150 ft. As 

shown in table 1 the 2 o (95 percent) arrival time error dispersion of nine of 
the eleven 4-D descents was within 19 seconds. Results of two descents were 
excluded because of defects in the experimental procedures; their arrival time 
errors were 25 and 41 seconds. It is expected that, in the future, the 2a 
dispersion can be reduced to about eight seconds by means of improvements in 
wind modeling, aircraft /engine modeling, navigation procedures and improved 
automation. These topics will be discussed more fully in the body of the 
report. Position and altitude errors, as shown in table 1, are well within 
the criteria for the flight tests. Airspeed errors at the end of descent were 
a little high, however, incorporation of the above four improvements is expected 
to maintain airspeed within ±5 knots at the end of descents. 

Preliminary analysis indicated that a linear wind model in which wind 
velocity is assumed to decrease linearly with altitude would be adequate for 
achieving end-of-descent arrival times within the ±30 second criteria for 
the worst expected wind conditions. Present thinking is that 4-D time control 
accuracy should be better than 30 seconds; this would indicate that a more 
exact wind model than the linear model is in order. 

Flight test wind velocity data were analyzed for two sets of tests 
wherein the aircraft performed descents over the same course between the 
altitudes of 33000 and 11000 ft. In the first set of descents the wind 
velocity profiles remained relatively unchanged for over five hours. In the 
second set .there was considerable variation of profile caused by frontal 


TABLE 1. - FLIGHT TEST RESULTS AND FUTURE PROJECTION 



1 

4-D Test 

Summary 

Future 

Error 


Goals 

(9 Descents) 

Projection 

End of Descent 

Mean 

Not specified 

-2s 

Os 

Arrival 

Time 

2 cr 

<30 s 

19s 

8s 

Airspeed 

±5 kts 

±1 5 kts 

+5 kts 

Position 

+1 n.mi. 

+0.2 n.mi.. 

±0.2 n.mi. 

Altitude 

+150 ft 

+50 ft 

+50 ft 
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activity in the area. These two sets of data provided useful insights in the 
comparison of a linear wind model and a more complex segmented model based on 
actual measurements. The segmented model improved wind model fidelity in the 
stable wind velocity situation but actually degraded fidelity in the situa- 
tion where the wind velocity profile was disturbed by the changing weather 
conditions. These two examples are not meant to provide a quantitative basis 
for wind modeling, however, the need to consider the effects of weather condi- 
tions on wind persistence is indicated, especially if actual wind measure- 
ments are the basis of wind profile modeling. 

Control laws for the 4-D descent were synthesized using a highly accurate 
non-linear model of the L-1011 on an interactive computer graphics terminal. 
Software development was carried out on an IBM 370 emulation of the ARMA 
flight management computer. Prior to flight the flight management computer 
was connected to a hot mockup in which the computer and its software were 
checked out with simulated aircraft system inputs. The flight tests were con- 
ducted without any serious complications, although, as mentioned above, two 
of the eleven descents were degraded by procedural problems. Remedies have 
been found that will be incorporated on downstream 4-D flight tests. 

The 4-D descent technique used in this study was a prototype technique 
that was tested to demonstrate the capability of achieving accurate arrival 
time control. The flight tests at the Dallas/Ft. Worth Regional Airport 
(DFW) were very instructive in learning how to refine the 4-D technique for 
the future so that it will fit smoothly into the air traffic control system 
and contribute to the solution of the air traffic congestion problem. A 
number of meetings were held with ATC personnel before and after the flight 
test demonstration at DFW; these sessions pointed out the need for the air- 
craft system developer to work toward the solution of flow control problems 
in cooperation with the ground controllers and the airlines in order to 
develop new procedures that are practical and that will be accepted by those 
who have to use them. 


INTRODUCTION 

This report documents the analysis performed and flight test results 
obtained by the Lockheed-Calif ornia Company under contract (NAS1-15546) to the 
National Aeronautics and Space Administration Langley Research Center at 
Hampton, Virginia. This effort was started in October 1978 with an analytical 
study of time-controlled (4-D) profile descent techniques and a review of wind 
modeling for the descent computation. The study was subsequently expanded to 
include flight tests on the Lockheed L-1011 which verified the feasibility of 
the 4-D high-profile descent. 

A motion film supplement is available on loan. A request card form and 
a description of the film are found at the back of this report. This film 
describes the flight to the Dallas/Ft. Worth Regional Airport which was the 
culmination of the L-1011 testing performed in the program. 
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LIST OF SYMBOLS 

Values are presented in both SI and U.S. Customary Units. Calculations 
were made in U.S. Customary Units. Altimeter readings are given in feet to 
correspond to the instrument reading. To convert from feet to meters multiply 
the value of feet by 0.3048. 

AC 


Alternating Current 

ACC 


Accessory 

ADI 


Attitude Direction Indicator 

AFCS 


Automatic Flight Control System 

A/L 


Autoland 

ARTCC 

Air Route Traffic Control Center 

ATC 


Air Traffic Control 

ATS 


Auto Throttle System 

B*D 


Beg inn ing-of -Descent 

CADC 


Central Air Data Computer 

CAS 


Calibrated Air Speed 

CDU 


Control and Display Unit 

CMD 


Command 

CPU 


Central Processor Unit 

CRT 


Cathode Ray Tube 

DA 


Drift Angle 

DC 


Direct Current 

DCAS 


Direct Computing Access System 

DFW 

4 


Dallas/Fort Worth Airport 

L . 





DME 

ECS 

E*D 

EPR 

ETA 

FAA 

FD 

FMC 

FMS 

GMT 

HSI 

HW 

IAS 

INAV 

INS 

I/O 

M 

NASA 

NAV 

NCU 

PMD 

RAD/RNAV 

RNAV 

S 

SAT 

STAR 

S/N 

TACAN 

TAS 

TCV 

TGT 

TH 

TKE 

TRACON 


Distance Measuring Equipment 
Environmental Control System 
End-of -Descent 
Engine Pressure Ratio 
Estimated Time of Arrival 
Federal Aviation Administration 
Flight Director 
Flight Management Computer 
Flight Management System 
Greenwich Mean Time 
Horizontal Situation Indicator 
Head Wind 

Indicated Air Speed 
Inertial Navigation 
Inertial Navigation System 
Input/Output 
Mach Number 

National Aeronautics and Space Administration 
Navigat ion 

Navigation Computer Unit 
Palmdale Airport (California) 

Radio/Area Navigation 

Area Navigation 

Laplace Operator 

Static Air Temperature 

Standard Terminal Arrival Route 

Serial Number 

Tactical Air Navigation 

True Air Speed 

Terminal-Configured Vehicle 

Turbine Gas Temperature 

True Heading 

Track Error 

Terminal Radar Approach Control 
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TW 

V 

Ve 

VGS 

VHF 

Vn 

V-NAV 

VOR 

VORTAC 

VT 

VW 

h 

h CMD 

h 

e 

r cmd 

r 

t 

a 

P 

P 

cr 


7 

Q 

0 CMD 

s 

Ah 

Ar 


At 

Atemp 

^sp 

^thr 


Tail Wind 

Velocity 

East Velocity 

Ground Speed 

Very High Frequency 

North Velocity 

Vertical Navigation 

VHF Omnidirectional Range 

Co-located VOR and TACAN 

True-air-speed 

Wind velocity 

Altitude 

Altitude Command 
Altitude Error 
Range Command 
Range 
Time 

Angle of Attack 
Regression Coefficient 
Correlation Coefficient 
Standard Deviation 
Flight Path Angle 
Pitch Angle 
Pitch Angle Command 
Summation 

Incremental Altitude 
Incremental Range 
Incremental Time 

Deviation from Standard Air Temperature 
Incremental Spoiler Command 
Incremental Throttle Command 


1.0 DISCUSSION OF TECHNICAL APPROACH 

This section describes the development of a time-controlled (4-D) high 
profile descent capability and its mechanization using existing L-1011 flight 
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management system (FMS) hardware and software. The production configuration 
FMS computes an optimum descent trajectory for the L-1011 and provides auto- 
matic control capability to guide the aircraft from cruise to an end-of-descent 
point (E*D) , arriving at a desired altitude and speed. To provide a 4— D 
capability a time constraint was added to the E*D arrival requirements by 
means of software changes to the existing computer program. 

1.1 The L-1011 Flight Management System 

1.1.1 Basic operation . - The FMS is an extension of the area navigation 
(or RNAV) capability originally certified with the aircraft in 1971. It per- 
forms the basic RNAV functions of waypoint navigation and coupled guidance 

as well as the automatic selection of VORTAC stations, tuning of the aircraft's 
VOR/DME receivers, and the mixing of inertial, radio, heading and air data 
sensor inputs to provide optimal navigation accuracy and reversionary mode 
operation in the event of degradation. The system is comprised of a computer, 
a CRT control and display unit (CDU) , and a CRT map display. 

The flight management capabilities, mainly related to automatic control 
of engine performance for all phases of flight, are as follows: 

• selection and performance of fixed or calculated optimum climb speed 
schedules 

• selection and performance of desired engine pressure ratio (EPR) with 
automatic or manual derating 

m calculation of optimum cruise conditions (i.e., altitude, speed, 
step-climb determinations) with automatic transitioning from climb 
to optimum or manually-specified cruise flight 

• calculation of the descent trajectory required for optimum or 
manually-specified descent speed schedules with automatic initiation 
and termination of descent at a prescribed end-of-descent point, at 
the desired altitude and speed 

• other capabilities such as estimated time enroute to waypoints, engine- 
out drift down modes, calculation of flap holding speeds, reversionary 
airport fuel and point-of-no-return calculations based on the effects 
of current and forecast winds which can be entered into the system by 
the crew. 

A description of the L-1011 FMS optimization algorithms is given in 
reference 1. 

1.1.2 System interface . - The typical aircraft configuration is a dual 
system installation with the electronic map accepting inputs from either 
system as desired. 

The functional interface of each system with the aircraft is summarized 
in figure 2. The major aircraft systems involved are the: 
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Figure 2. - Flight management system block diagram 


















• inertial and radio navigation sensors 

• heading reference system 

• air data systems 

• automatic flight control and flight director system pitch and roll 
channels 

• automatic throttle 

• related flight instruments. 

The detailed interface, shown in figure 3, illustrates that considerable 
redundancy exists for the air data and navigation functions. Single system 
failures are automatically dealt with and the most optimum operational 
mechanization is configured for the sensors available at any given time. The 
available navigation modes are: 

• inertial /radio mix with up to three inertial system and two each 
VOR/DME receivers 

• inertial only (one to three inertial systems) 

• radio only with sub modes of 

• DME/DME 

* VOR/DME 

• heading/air data 


1.2 Descent Profile Modeling 

A precomputed descent profile is the basis of the 4-D control system 
synthesized in this study. This 3-D descent profile is computed by the FMS 
prior to descent using latest aircraft weight, wind and temperature data. It 
is based on descending at near idle power according to one of several available 
fixed Mach/IAS speed schedules. In the 4-D system designed for this program 
feedback loops are added to the 3-D system so that variations in range and 
altitude from the precomputed trajectory are corrected by activating the stabi- 
lizer, throttles, or speed brakes. The 3-D profile is stored in the FMS com- 
puter as a function of time. At each point in time during the descent the 
spatial position of the aircraft is compared with the stored trajectory to 
evaluate the altitude and range errors. The paragraphs which follow describe 
the steps taken in the development of this capability. 


1.2.1 Descent model verification. - At the outset of the study the 
fidelity of the 3-D descent computation made by the actual airborne FMS was 
compared with actual flight test data of the L-1011 in descent. Thirty-six 
descents were simulated in the Palmdale hot mock-up lab facility using the 
latest production FMS hardware and software. Altitude vs. range data were 
taken for various combinations of initial conditions such as descent speed 
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Figure 3. 
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interface . 
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schedule, aircraft weight, temperature, and end-of-descent point altitude. 
These data were then entered into the interactive graphic computer simulation 
program at Burbank and compared with the latest L-1011 flight test data. 
Refinements to the FMS 3-D algorithms were then specified for the airborne 
computer program. 


1.2.2 Computation of the 3-D descent profile . - The FMS 3-D descent 
profile is generated in six segments, back-computed from the specified end- 
of-descent point (E*D) , as shown in figure 4. The profile is established by 
sequential build-up of incremental ranges ( Ar^) and altitudes ( Ah^) for each 
of the six segments. Ah is set to a constant value of 500 feet except for 
the level flight segments, and for segments transitioning to or from level 
flight. The corresponding Ar's are then computed as a function of Ah using 
pre-stored polynomials determined from L-1011 aerodynamic data. The point at 
which the back-computed profile intersects cruise altitude is the beginning- 
of-descent point, or B*D. The calculated position of B*D is then automatically 
entered as a waypoint by the FMS computer program into the flight plan previ- 
ously selected by the crew. The descent profile and the position of B*D are 
continually updated in flight to reflect current predictions of the descent 
entry initial conditions. The definition of B*D as a waypoint allows the 
program to calculate and display other desirable flight plan-related informa- 
tion, such as estimated t ime-of-arrival (ETA), course, distance, and time 
enroute to B*D. 


1.2.3 Computation of time-in-descent . - In order to 
of arrival at E*D, the total time-in-descent is computed. 
E*D can be computed by simply adding the calculated total 
the estimated time of arrival at B*D: 


predict the time 
Then the ETA at 
time-in-descent to 


E*D ETA = B*D ETA + calculated total time-in-descent 


Total time-in-descent is determined by first calculating an incremental 
time (At-^) required to descend through each Ah-^ /Ar^ segment of the trajectory; 
the total time-in-descent is then the summation of these incremental times or 


Time-in-descent = SAt^ 


The basic relationships used in the At^ computation are as depicted in 
figure 5. 

With wind velocity (VW-^) , true-air speed (VT-^) , Ah^ and Ar^ given at the 
i-th altitude level, ground speed (VGS-^) and At^ are then computed: 
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VGS . 
1 


At 


VT . 
1 


cos y . + vw. 

X X 

Ar . 


x 


1/2 (VGS. + VGS. , _ . 

x x+1 ; 

2 A r . 
x 


(VT. + VT . ... ) Cos y . + vw. + vw 
x x+1 x x 


i+1 


where CosY. is approximated by the earth referenced path angle 


Cos Y . ss Cos "Y 
x e 


A ri 


V Ah i +Ar J 


y : 

e 


This approximation, used because the wind-referenced path angle 7^ 
was not calculated by the existing FMS software, results in an open loop error 
of 2 seconds or less, and therefore had negligible effect on the closed-loop 
performance of the prototype system. 

The time-in-descent algorithms developed for all six segments of the 
descent profile are included in Appendix A. 

1.2.4 Descent profile sensitivities and error sources. - The effects of 
speed schedule, aircraft weight, temperature and E*D altitude on descent 
trajectory range and time- in-descent were analyzed. The results of this 
analysis are summarized in figures 6 and 7 and presented in more detail in 
Appendix B. 

Error sources which can significantly affect the quality of the descent 
trajectory and time-in-descent models were determined to be: 

• aircraft performance modeling 

• engine modeling 

• navigation sensors 

• modeling of winds aloft and along the descent 

The aerodynamic performance model used for this program was the latest 
available for production aircraft. Figure 8 illustrates the operational quality 
of the 3-D profile using this model. No attempt was made to tailor the trajec- 
tory computations specifically for the flight test aircraft, S/N 1001, because 
it is continually being used for a wide variety of flight test activities in- 
volving engine testing and the development of other advanced systems (e.g.. 
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Distance from B*D (n.mi.) 

Figure 8. - Descent profile. Test No. 216, 8/6/79. 

active controls) . It is not uncommon to fly with a "mix" of various engine 
types and various control surface configurations, which are subject to change 
on short notice. Producing a more flexible model, that could be modified as 
the aircraft configuration was changed, was not possible in the short span of 
the study. It was noted that the engine mixes, of which there were several 
during the program, and the presence of wing tip extensions, that were not 
represented in the production modeling, made appreciable differences in actual 
descent trajectories. The closed loop nature of the 4-D control system sup- 
pressed errors from these sources. For future flight testing it is, however, 
recommended that more exact modeling be used. 

The FMS inertial/radio navigation solution with DME/DME updating which 
offers position accuracies of 0.5 n.mi. or better was considered to be more 
than adequate for this application. However, as discussed in a later section 
of this report, further consideration concerning the logic utilized in blend- 
ing the navigation sensor "mixes" appears necessary to preclude the position 
updating processes from introducing apparent time errors during 4-D flight. 

Wind has the potential of causing the largest 4-D descent performance 
errors. Therefore it is necessary to accurately measure the wind existing in 
cruise just prior to descent, and to construct a realistic model for the winds 
existing at lower altitudes, so that the descent trajectory can be accurately 
computed. The measurement of winds in cruise presented no problem; true 
airspeed, ground speed, ground track and heading were used for this calculation. 
The descent wind model used in the FMS during flight testing was a linear decay 
of the wind velocity at cruise to zero at sea level. 
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1.3 Wind Study 


The accurate modeling of winds along the descent trajectory is a critical 
factor in the computation of the descent profile and in the calculation of the 
t ime- in-descent . The following paragraphs describe a study of existing wind 
data and the basis for selecting the linear wind model that was used in the 
flight test experiments. 


1.3.1 Statistical wind data . - In earlier NASA studies wind data (see 
references 2 and 3) were compiled for seven geographical locations of the 
world. The data, collected over a period of several years at each location, 
were presented as monthly tabulations of the zonal (east-west) and meridional 
(north-south) mean wind components at integral kilometers of altitude. Stan- 
dard deviations from the mean were also presented for each altitude, along 
with the correlation coefficients relating the winds at one level upon another. 
Table 2 is an example of this data for Tripoli, Libya, during the month of 
February . 

In the present study the monthly data from each location, up to the certi- 
fied operational ceiling altitude of the L-1011 (42000 ft), were entered into 
an IBM 370 computer at Burbank for additional analysis. This effort is de- 
scribed in paragraphs 1.3.2 through 1.3.4. 


1.3.2 Linear wind model development . - Reference 4, reporting the 
results of a study concerning the correlation of inter level wind velocities, 
states that a linear relationship exists between winds at cruise and the winds 
along the lower altitudes. This relationship is illustrated by the scatter 
diagrams presented in figure 9 where linear regression techniques were used to 
provide a best fit to the observed interlevel wind velocity data. The slopes 
and intercepts of these regression lines can be determined using conditional 
probability density functions; the resulting linear equations can then be used 
to compute the expected values of wind velocities of all lower altitudes when 
the winds at cruise are known. 


V 

i 


= V. + 

l 




(V. 

1 


V.) 

3 


where 


(3. . = P. . 
iJ ij 


cr. 
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TABLE 2 


NASA WIND DATA 










































Wind speed at lower altitude (Knots) 






is the regression coefficient; i =0,1,2--- (j-1) are the lower altitude 
levels at which the expected values of wind velocities are to be computed; 
j is the altitude level at which the known wind velocity is given; is the 
computed wind velocity at level i; Vj is the known wind velocity measured 
at level j; V-j_ and Vj are the mean wind velocities; P-jj is the interlevel 
correlation coefficient for level j on level i; c t± and cr-j are the standard 
deviations. The 95-percentile confidence envelope for the conditional wind 
profile is given by: 


V.' = V. ±2it. 1 

it i 



1/2 


A computer program was written for the IBM 370 to investigate the 
possible application of these techniques in developing a suitable wind model 
for the descent trajectory computation. Wind profiles based on the linear 
regression techniques were computed and plotted for each of these locations. 

Figure 10 shows the February zonal mean wind profile together with the 
+2cr raw data envelope and the +2cr conditional profile envelope. Figure 11 
shows the February zonal mean profile and the conditional profiles resulting 
from +lcr and +2(r winds existing at a descent entry altitude of 42000 feet. 

A review of the graphic results concluded that, for the flight regime 
of the L-1011 aircraft, wind velocities at cruise altitudes decay in a fairly 
linear fashion to zero knots at zero feet altitude. 


1.3.3 Simulation wind model development . - For the seven geographical 
locations studied, winds were observed to be predominantly east-west; the 
worst case wind (i.e. the monthly wind velocity with the highest standard 
deviation) was seen to occur at Tripoli, Libya in the month of February. For 
this reason, the February zonal Tripoli wind profiles were chosen for use as 
a worst case statistical representation for simulation studies in the course 
of development of the 4-D descent control algorithm. To add realisn to the 
descent simulations using the Tripoli model, a random component, taken from 
Palmdale flight test data, was added. The Palmdale wind, with a standard 
deviation of 7.6 kts, was added to the Tripoli wind as shown in figure 12. 

The resulting wind was then applied as a headwind or a tailwind during the 
analyses . 


1.3.4 The effect of wind estimation error on open loop performance . - 
The following descent performance criteria were established at the beginning 
of the program: the E*D position error goal was ±1 n.mi. and E*D time-of- 

arrival error was to be kept within ±30 seconds (2 cr) . The Tripoli 2<r raw 
wind profile and the mean wind profile, shown in figure 10, were applied as 
headwinds to the descent profile algorithm described in section 1.2. The 
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Figure 12. 


Simulation wind models. 


differences in the descent time and descent range values resulting from the 
application of these two wind profiles represent the open-loop 3-D descent 
errors that can result from a worst-case unmodeled wind. The results of this 
study showed that E*D could be overflown within about ±40 seconds of the re- 
quired time > however, E*D altitude would be captured with range errors of 
approximately ±20 n.mi. This is presented by figure 13 which illustrates the 
inadequacy of an open loop solution for achieving the desired results. 


1.4 4-D Control Law Development 


1.4.1 Altitude and rang e errors as feedback variables . - With the need for 
improved 4-D control clearly justified by simulation results, several candidate 
schemes were identified and traded-off. The technique eventually chosen was to: 

• Precompute a fixed 3-D trajectory prior to descent 

• Feedback trajectory altitude errors to the spoilers or throttle 

• Feedback trajectory range errors through the aircraft pitch autopilot 

Figure 14 depicts the feedback variables of altitude error and range error 
that are used for control. 

Control of altitude errors was not automated because of the extensive 
hardware modifications that would be required for the aircraft : s flight con- 
trol systems. Instead, advisories of altitude deviations from the 3-D profile 
were displayed to the flight crew for manual deployment of spoilers (to 
increase the rate of descent for the "too high" case) or manual application of 
thrust (to slow the descent rate for the "too low" case). The range error 
feedback was completely automatic through the pitch autopilot. 

It was decided to apply control only during the airspeed hold portion of 
the descent (below 29,000 ft altitude) for two reasons: simulation studies 

showed that altitude and range (time) errors that might accumulate early in 
the Mach-hold portion of descent were not significant and could be compensated 
for by applying control in the airspeed region; the Mach-hold control laws and 
logic in the existing FMS are fairly complex and not easily modified by 
additional software. 

One descent speed schedule was mechanized for the prototype system. This 
schedule consisted of a 0.82 Mach-hold descent until indicated airspeed (IAS) 
increases in value to 320 kts (which occurs at approximately 29000 ft); after 
this, the system holds 320 kts IAS until level-off at 10000 ft or at E*D 
altitude, whichever comes first. A 0.82 mach/320 kts descent was selected 
primarily because a ±40 kts airspeed variation would be available for closed 
loop control of range errors. It also represented fairly closely a minimum 
cost speed schedule. 
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Figure 14. - Feedback variables. 


1.4. 1.1 Altitude error feedback . - Figure 15 shows the manual control 
scheme utilized for making good the 3-D or spatial profile. The aircraft’s 
actual altitude is compared to the altitude required by the descent profile 
to derive altitude error. For positive errors, i.e., when the aircraft is 
above the desired spatial profile, spoilers are deployed to increase the air- 
craft's descent rate; conversely, thrust is applied to slow the aircraft’s 
descent when the aircraft is below the profile. A "bang-bang" control strat- 
egy was adopted, wherein errors were allowed to build up to ±300 ft and speci 
fied spoiler/EPR schedules were followed to reduce the error to ±150 ft. Thi 
was done because a proportional control strategy: 

• would have been subject to interpretation by the crew and resulting 
performance would have been more difficult to evaluate by simulation 
analysis 

• would have caused more throttle cycling than the bang-bang approach 

The EPR schedule shown in figure 15 was calculated for an L-1011 with a 
full complement of RB 211-524B engines. 
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Figure 15. - Altitude error feedback. 
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1.4. 1.2 Range error feedback . - The time control loop shown in figure 16 
is switched in as the aircraft transitions from mach hold to airspeed hold in 
the descent, for a speed schedule of 0.82 mach/320 kts, this occurs nominally 
at 29000 ft. Control remains active to 10000 ft or until the aircraft cap- 
tures E*D altitude, whichever comes first. Gains were chosen to allow the 
airspeed to vary slowly (0.2 kts/sec) between 280 and 360 kts (320 ±40 kts) 
during airspeed "hold” to make good the E*D arrival time requirement. 


1.4.2 Descent simulation errors . - Several simulation runs were made to 
examine the performance of the altitude and range control strategies discussed 
above. The wind profiles shown in figure 12 were used to introduce unmodeled 
wind errors during descent; a three degree-of-f reedom dynamic aircraft model 
was used to evaluate the effectiveness of the 4-D control laws. Table 3 tabu- 
lates the errors in time, range, and E*D altitude/speed resulting from these 
simulations. Time response plots are presented in Appendix C. 


TABLE 3. - DESCENT SIMULATION ERRORS 


Errors at 
End of 
Descent 

4-D Control Off 

4-D Control On 

Time 

(s) 

Altitude 

(ft) 

IAS 

(kts) 

Range 

(n.mi.) 

Time 

(s) 

Altitude 

(ft) 

IAS 

(kts) 

Range 

(n.mi.) 

Palmdale Wind 

-4.3 

97.6 

3.4 

0.33 

3.5 

16.7 

3.3 

-0.36 

Palmdale Wind 
+ 

Head Wind 

59.1 

0 

2.2 

-4.2 

22.0 

11.9 

3.4 

-1.6 

Palmdale Wind 
+ 

Tail Wind 

-64.3 

1000 

11.3 

5.4 

-15.7 

141 

5.2 

1.3 


Note: Time, altitude, and IAS errors were measured when the aircraft position was at the metering fix. 


Range error was measured at the required-time-of-arrival at the metering fix. 


1.5 FMS Software Development 

The actual preparation of software incorporating the 4-D descent capa- 
bility into the L-1011 Flight Management System was funded internally by 
Lockheed because of the proprietary nature of the existing FMS algorithms and 
control laws. This report therefore contains no detailed software coding 
or other related information in order to protect the rights of Lockheed and 
ARMA, Division of AMBAC, the system's manufacturer. The software develop- 
mental approach is discussed, however. 
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Figure 16. - Range error feedback. 





1.5.1 Computer organization . - The L-1011 Flight Management System 
computer used for this work was an ARMA model 1808 with 32K words (18 bit/word) 
of core memory; the CPU has 56 different instructions. The functional charac- 
teristics of the ARMA 1808 are: 

• .909 MHz clock rate 

• 700 psec memory access time 

• fixed point 2's complement arithmetic 

• execution time: 6.6 p sec add/subtract 

26.4 p sec multiply /divide 

• programming registers: 1 upper accumulator (A) 

1 lower accumulator (B) 

3 index (XI, X2 , X3) 

1 program counter 

• address modes: direct 

indirect 

indexed 

• A/D and D/A conversion 

• Synchro I/O 

• A/C DC I/O 

• VOR/DME auto tune 

• Discrete I/O 

• CDU keyboard interface 

• Serial digital I/O channels 

• Direct memory access (up to 151,000 words/sec) 

• Priority interrupt 

• 4.5 msec real time counter 


1.5.2 Emulator development . - The 4-D descent software was written as 
machine language patch code to the basic Saudia Airlines FMS operational pro- 
gram assembly. Coding relating to the FMS engine-out mode was removed and the 4-D 
patch code was developed in this section of core, thereby preserving the rest of 
the FMS functional capability. 

A FORTRAN IV emulation, reproducing the instruction set of the ARMA 1808 
on a bit-by-bit basis, was prepared and used for the initial development of 
the 4-D software code. In this way, code could be written in the same format 
required by the airborne computer and tested statically in an interactive en- 
vironment. The code could be stepped through, one instruction at a time, or 
the program executed with various check points in the flow; the results of 
any operation could be displayed in engineering units rather than fixed-point 
binary, and logic and branching processes easily reviewed. 
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1.5.3 Data bank preparation . - The Flight Management System data bank 
is a section of memory in the computer which is reserved for storing necessary 
information primarily related to the navigation function of the system. Exam- 
ples of information stored in the data bank include: 

• waypoint names and locations 

• VORTAC data 

• route structures 

Data for the Dallas /Ft. Worth area were collected, and a punched computer 
tape made in preparation for the demonstration flights that were performed on 
01 August, 1979. Appendix D presents this information in detail. 


1.5.4 Hot mock-up checkout. - After the software changes to the opera- 
tional program had been checked out on the emulator and the data bank informa- 
tion prepared, the machine code patches were transferred to the airborne com- 
puter and system level testing conducted on the Palmdale Hot Mock-up facility; 
this facility allows the airborne system and its interface to be exercised in 
a dynamic simulation of the aircraft. In this way, the total FMS software 
package and its interaction with the 4-D patch code was evaluated prior to in- 
stallation and flight testing in the Lockheed L-1011 test aircraft. 

The overall process of software development is illustrated in figure 17. 



Figure 17. - Software development. 
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2.0 FLIGHT TESTING 


Flight testing of the new '4-D descent capability was conducted on L-1011 
S /N 1001, shown in figure 18. A total of eleven descents were performed - 
eight in the local area around the Flight Test Center at Palmdale, California 
and off the Pacific coast, and three in the arrival air space of the Dallas/ 
Ft. Worth Regional Airport (DFW) . Approximately sixteen flight hours were 
used in the flight test evaluation of the 4-D descent. 


2.1 Local Area Flight Testing 

The flight test effort began on June 25, 1979, with two straight-in des- 
cents performed between 31000 and 11000 ft to minimize any possible errors 
that might be introduced by turns in the descent and by extended operation in 
the mach-hold region. 

Two more straight-in descents were performed on the 26th and the 29th - 
the first between 37000 and 11000 ft to test the exposure to uncontrolled 
flight in the mach region, the second between 31000 and 9000 ft to test 
the exposure to uncontrolled flight at altitudes below to 10000 ft. 

Four descents took place on July 30 in final preparation for the Dallas 
demonstration flights. The descents were conducted off the Pacific coast 
over an altitude range and a course which simulated the approach geometry 
of the Acton Standard Terminal Arrival Route (STAR) at DFW. 


2.2 Dallas/Ft Worth Demonstration Flights 

The last three descents in the flight test program were performed as 
part of a demonstration flight into DFW on August 1, 1979. Figure 19 shows 
the flow of air traffic to DFW for final approach and landing to the south. 
Four "corner post" VORTACS are used as the metering fix gateways defining the 
approximate boundaries between the Ft. Worth Center and Dallas approach con- 
trol. All three descents were performed from the southwest, via the Acton 
Six STAR. The waypoint Flato was chosen as E*D with an altitude and speed 
requirement of 11000 ft and 250 kts . 

Figure 20 illustrates the westerly approach (Wink transition) to the 
Acton STAR that was flown for the first demonstration. For this first descent, 
B*D was calculated co lie between Tuscola and Gramy at a point along the 
flight plan 92.4 n.mi. before Flato. Although not a primary objective of the 
demonstration, Wink, Texas, was overflown within 30 seconds of the time re- 
quested by Ft. Worth Center prior to the aircraft's departure from Palmdale, 
California. 

The L-1011 FMS had no capability for performing controlled (closed loop) 
4-D cruise flight, however, full use was made of the system’s preflight 
planning capability to achieve a very credible approximation to 4-D flight 
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Figure 18. - Lockheed Advanced TriStar, L-1011 S/N 1001 
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Figure 20. - Acton arrival Wink transition. 



management. Prior to takeoff at Palmdale, Lockheed Jet Plan wind and 
temperature data were entered into the FMS for the cruise portion of flight 
to obtain an accurate projection of time enroute to Wink. The take-off time 
was then backcomputed to meet the 13:00 Wink arrival time requirement. The 
L-1011 took off at this time, 09: 17 local, and flew to Wink under automatic 
FMS control. 

Figure 21 is a schematic representation of the demonstration flight 
scenario. The first high profile descent, using the inertial navigation system 
as a reference, was performed at the end of the flight which originated in 
Palmdale, California. The aircraft, still under automatic control, coupled 
to the ILS for final approach and landing. After landing, observers from 
NASA, FAA,and the airlines were boarded for two more descent demonstrations. 

An inertial/radio mix was used for these descents for improved accuracy. After 
the third descent, the aircraft again made a fully automatic approach and 
landing. The day after the flight, a debriefing workshop was held in which 
representatives from NASA, the FAA, Delta Airlines, and Lockheed discussed 
the results of the demonstrations. 

Table 4 lists in detail the initial conditions for each of the eleven descents, 
performed in the course of the flight test effort. 


TABLE 4. - FLIGHT TEST INITIAL CONDITIONS 






R9 






r 

~ 



Altitude 

Wind at 

Mmw 

SAT 

Descent 

Descent 

Nav 

NCU 

Descent 

Date 

• B*D/E*D (kft) 

B*D(kts) 

in 

(°C> 

range (n.mi.) 

Time (s) 

mode 

S/N 

■1 

6/25/79 

31/11 

46 HW 

138 300 

-37 

68.7 

568 

INERTIAL/ 

251 


6/25/79 

31/11 

38 TW 

126 300 

-39 

74.6 

571 

RADIO 

1 




6/26/79 

37/11 

56 TW 

136 000 

-54 

98.1 

761 





111 

6/29/79 

31/9 

33 TW 

143 800 

-38 

85.7 

731 

INERTIAL 



1 

7/30/79 

31/11 

19 HW 

173 275 

-32 

81.5 

633 



21 

03 

mm 

7/30/79 

31/11 

26 HW 

160 800 

-32 

76.5 

604 





1 

7/30/79 

31/11 

31 HW 

150 595 

-32 

73.4 

581 





8 

7/30/79 

31/11 

24 HW 

133 980 

-32 

71.5 

554 





9 

8/1/79 

. I 

27 TW 

171 000 

-36 

92.4 

693 

INER 

TIAL/ 



10 

8/1/79 


28 TW 

161 000 

-36 

88.1 

686 

RADIO 



11 

8/1/79 

33/11 

29 TW 

150 000 

-36 

87.0 

658 

1 




Note: Speed schedule for all descents - 0.820/320/250 
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2.3 Flight Test Results 

As shown in Table 5 the most important flight test goals, arrival time, 
airspeed, position, and altitude accuracies, were mainly achieved; arrival 
time accuracy is of greatest interest in a 4-D experiment. It was an objective 
to arrive at the end-of-descent point within ±30 seconds of the predicted time; 
this goal was achieved with a 2 tr error of 19 seconds about a mean of -2 seconds. 
Figure 22 is a histogram of all eleven descents reflecting the even distribu- 
tion of E*D errors based on the last estimate made by the FMS computer before 
descent. The values in Table 5 were obtained from nine of the eleven descents; 
two were eliminated because of test procedure problems. The position and 
altitude errors of Table 5 are also within the pre-established test goals. 
Airspeed error was high, i.e., within ±15 kts instead of ±5. This amount of 
error is not substantial and will be reduced with system improvements dis- 
cussed in paragraph 2.5. 

The terminal errors observed at E*D for each of the eleven descents are 
presented in table 6. Time errors presented under the column heading "Initial 
Estimate" were computed by the FMS in cruise, prior to B*D, as described in 
paragraph 1.2.3. The arrival time estimate for E*D given in this column 
would vary with changing winds. There was concern therefore that the predicted 
time of arrival given to ATC by the crew about 5 minutes before descent might 
not be accurate, however, this concern was unfounded. Mean and 2a values in 
table 7 for the initial and final estimates agree within one second. 

The method used to measure terminal errors at the end-of-descent was to 
analyze recorded time, altitude, and airspeed for the interval when E*D was 
overflown. System position error was not measured against a ground-based ref- 
erence for the first four descents. DME range was used as a position refer- 
ence for descents five through eight; approach control radar data provided the 
position reference for the last three descents, which were performed at Dallas. 

The pre-selected E*D altitude is captured by the autopilot prior to E*D 
as the aircraft begins to decelerate to the required final speed; altitude 
errors at E*D are therefore determined by the accuracy of the autopilot's 
altitude hold mode and are typically less than ±50 feet. E*D speed is con- 
trolled by the autothrottle system using commands from the FMS. Unfortunately, 
a wiring problem which caused the autothrottle to disconnect at E*D altitude 
capture, existed on the aircraft during the period the first eight descents 
were performed; as a result, the autothrottle could not automatically capture 
and hold the required speed of 250 kts at E^D. The problem was corrected, 
however, in time for the last three descents performed at Dallas. 

Two descents were degraded by procedural problems. The third descent, 
performed on June 26, 1979, was initiated over the Pacific Ocean at a large 
distance from the VORTAC station being used as a navigation reference. As the 
descent progressed toward the mainland, radio position updating caused the 
position of the descending aircraft to change and the spatial path of the 
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TABLE 5. - FLIGHT TEST RESULTS 


Error 

4-D Test 
Goals 

Summary 
(9 Descents) 

End of Descent 

Mean 

Not specified 

2s Early 

Arrival 




Time 

2 or 

<30 s 

19s 

Airspeed 

±5 kts 

±1 5 kts 

Position 

+1 n.mi. 

+0.2 n.mi. 

Altitude 

±150 ft 

±50 ft 


Number of 
occurrences 



Seconds early M Seconds late 


Figure 22. - E*D time error histogram 
(final estimate). 

descent to move towards the updated E*D position. It is estimated that 20 
of the 25 seconds of time error observed could be attributed to this cause. 

The largest time error recorded during the flight test program occurred on 
July 30,1979, during the eighth descent. Data analysis indicates that 
incorrect data were entered into the system prior to descent. Aircraft weight 
was manually input to the computer as an expedient during this test program; 
it is planned to use an automatic weight computation for future 4-D testing 
as is the practice on production FMS systems. 

Figure 23 depicts the aircraft ground track for the first two descents 
into the Dallas terminal area on August 1, 1979, as recorded by approach 
control radar. The first descent was performed using an inertial system as 
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TABLE 6. - FLIGHT TEST TERMINAL ERRORS AT E*D 


Descent 

Time 

(Initial Estimate) 
(s) 

Time 

(Pinal estimate) 
(s) 

Altitude 

(ft) 

Speed 

(kts) 

Position 

Ref 

Error 

(n.mij 

1 

No estimate made 

4 Early 


60 Fast> 


RNAV 

Not n 

measured 

2 



8 Late 


■ 

0 






3 



25 Early (B) 


n 

0 






4 



14 Early 



0 

\/ A\ 



1 


5 

3 Early 


2 Late 


■ 

45 Fast 


DM 

E 

<2 

6 

13 Early 


16 Early 


■ 

50 Fast 






7 

3 Early 


6 Early 


■ 

20 Fast 






8 

66 Late (C) 

41 Late (C) 


■ 

0 






9 

14 Late 


12 Late 


■ 

10 Fast 


RA 

DAR 

2.5 

10 

3 Early 


6 Early 


■ 

15 Fast 




<0.2 

11 

1 Late 


5 Late 



8 Fast 





<0.2 


Notes: 


(A) Autothrottle did not automatically capture and hold 250 kts at E * D for descents 1 through 8 

(B) Most of this error caused byradio position updating during descent 

(C) Caused by incorrect input data to the system. 


TABLE 7. - ARRIVAL TIME ERROR COMPARISON 



Initial Estimate 

Final Estimate 


Error (6 Descents) 

Error (9 Descents) 

Mean 

Is Early 

2s Early 

2o* 

18s 

19s 


41 




























the navigation reference; radar data shows the aircraft's position to have a 
southerly error of 2.5 n.mi. The inertial system had not been updated since 
its preflight check nearly three hours before the end of descent. The second 
and third descents of the day were performed using radio reference navigation 
for improved accuracy; radar data shows the errors observed during these 
descents to be less than 0.2 n.mi. 


2.3.1 Wind modeling. — On two sets of flights, wind profile data were 
collected for several descents performed over the same course. Figure 24 shows 
wind profile data for the 7/30/79 tests in which four descents were made off 
the California coast starting at Saint Nicholas Island and ending over the 
Catalina Channel. The aircraft encountered descent headwinds which varied 
little in the 5 hours 29 minutes between the first and last descents. Fig- 
ure 25 illustrates data from the three descents performed 8/1/79 on the Acton 
STAR to DFW. These winds were predominately tailwinds and are seen to be less 
correlated than those of figure 24; the diminished correlation was caused by 
frontal activity and weather cells in the vicinity of the approach path. These 
two sets oi data afforded an opportunity to compare the system linear wind 
model with a segmented model. (the linear model was the one used for the 
flight test experiments as described in paragraph 1.3.2). The segmented 
model is one in which actual wind measurements made at three thousand foot 
altitude increments are connected by straight line segments which could then 
be used as a model for subsequent descents. These measurements could be 
collected at an airport from meterological instruments or from descending air- 
craft and radio-linked up as segments to other approaching aircraft. 

Figure 26 shows a segment fit to the winds encountered during the first 
descent on 7/30/79, similarly, figure 27 shows the segment fit to the first 
descent of 8/1/79. These segments were then plotted against the actual winds 
encountered during the last descent of each day. 

These plots were reconstructed with wind velocity on the ordinate and 
time on the abscissa as the right hand graphs of figures 28 and 29, also 
shown on figures 28 and 29, in the left hand graphs, are the linear models 
programmed for the last descent of each day superimposed on the actual wind 
profiles. (Note that the linear models have breaks to accommodate the course 
change made during the descents.) The areas between the curves on each graph 
were integrated to obtain a measure of the trajectory error, expressed in 
nautical miles, that would accrue during these descents for the respective 
models; these errors are summarized in table 8. 

As expected for the well-correlated 7/30/79 data, the segmented model 
has less absolute error than the linear model (0.16 n.mi. and 1.13 n.mi. re- 
spectively). For the descents performed on 8/1/79, the segmented model did 
not improve matters, in fact its error was worse (-1.9 n.mi. compared to 
-0.13 n.mi.). These two data sets are not meant to provide quantitative guide- 
lines for a phenomena that is best described by statistical means; they do, 
however, suggest some possibilities for future research. The persistence of 
wind velocity profiles in unstable weather conditions is such an area. 
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'lgure 26. - Segment fit to 7/30/79 

wind profile, first descent. 


Figure 27. - Segment fit to 8/1/79 

wind profile, first descent, 
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Figure 29. - Comparison of linear and segmented wind models, 8/1/79 flight tests, L-1011 S/N1001. 





TABLE 8. - WIND MODEL ERROR 



Flight test date 

7/30/79 

8/1/79 

Linear model 

1.13 n.mi. 

— .13 n.mi. 

Segmented model 

.16 n.mi. 

—1.9 n.mi. 

Age of segmented model 

5 hrs, 29 min 

3 hrs, 20 min 

-Denotes tail wind error 



+ Denotes head wind error 




Wind persistence in stable conditions has been researched and it has 
been shown to be quite good for periods up to six hours. (See reference 5.) 
The 7/30/79 data in figure 24 are consistent with this finding. 


2.3.2 Effects of wind modeling errors on 4-D descent . - Unmodeled wind 
errors are corrected by aircraft speed changes. The amount of speed change 
possible varies widely as a function of the descent speed schedule flown. 
Changes in E*D arrival time are also effected by varying the descent speed. 

A practical 4-D descent system must have the flexibility to vary its E*D 
arrival time somewhat to properly interface with the ground controller. This 
flexibility to vary arrival time without a change in flight path is eroded 
by unmodeled wind errors. The 4-D systems of the future therefore will require 
accurate wind models that will not cause a serious degradation of arrival time 
flexibility . 


2.4 Air Traffic Control Integration 

The development of 4-D techniques must be pursued with the goal of even- 
tually integrating them into the national air traffic control system. Based 
on the results of this study and the NASA TCV flight tests at Denver, it 
appears as if that goal is possible without major changes in ATC methods. 

(See references 6, 7, and 8.) DEW and Denver both make use of automated time- 
based metering, a method of local flow control that will be extended to other 
busy airports in the next few years. The integration of 4-D will be discussed 
from the standpoint of interfacing with such an automated metering system. 
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2.4.1 Metering and spacing . - Time-based metering uses computer-aided 
techniques to control the flow of aircraft into the terminal area at a rate 
that does not exceed the acceptance rate of the airport. Metering is the 
responsibility of the Air Route Traffic Control Center (ARTCC) . As shown 
in figure 30, spacing of aircraft in the vicinity of the airport is the respon- 
sibility of Terminal Radar Approach Control (or TRACON) . The metering program 
computes the arrival time of all aircraft approach the airport when they are 
about 25 minutes away. Arriving aircraft and their associated arrival times 
are then displayed on a CRT scope in their order of arrival for the ATC flow 
controller to determine if the rate of arrival will exceed the acceptance 
rate of the airport, if so, the arrival times are then delayed to match the 
current rate of acceptance. The various sector controllers are each given 
the metering fix arrival times assigned by the flow controller; they then 
radio speed or path changes to incoming aircraft to meet the assigned times 
for the aircrafts in their respective sectors. 


An aircraft with 4-D flight management capability can fit neatly into this 
system. To do so it is necessary for the ground controller to radio the 
assigned metering fix time to the aircraft so that the 4-D descent can be pro- 
grammed to cross the metering fix at exactly that time. This can be done now 
without any change in the metering program. The FMS computation would have to 
have sufficient flexibility to handle an assigned time that may not correspond 
exactly with the aircraft's ETA at the metering fix. At least thirty seconds 
of flexibility would be required to accommodate the fact that the metering fix 
times are given in integral minutes. Additional flexibility will be necessary 
to accommodate delays caused by weather or traffic peaks. 

Substantial reductions in ground controller workload are possible with 
4-D techniques; additionally, 4-D can do the job better. A controller can 
bring aircraft to the metering fix within one to two minutes of the assigned 
time whereas the 4-D equipped aircraft can get there within seconds of the 
assigned metering fix time. 

Extending 4-D back to takeoff has the potential of greatly reducing 
congestion. As was demonstrated in the DFW flight test, the L-1011 flew over 
Wink, Texas, within 30 seconds of the time assigned by ATC before takeoff. 

Often the same scheduled flights are delayed at the same time every day when 
approaching busy hub airports. If all scheduled aircraft were assigned arrival 
tines and used 4-D now, it is apparent that the delays which occur at peak 
periods could be reduced or eliminated because there wouldn't be bunching up 
at the metering fixes. 
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Figure 30, - Metering and spacing. 



The FAA is already using the Fuel Advisory Delay System wherein aircraft 
scheduled into a busy airport, such as O’Hare Airport in Chicago are advised 
to delay takeoff if their normal arrival occurs at a time when the airport 
acceptance rate is projected to be exceeded. Total 4-D flight planning would 
integrate well with this system also. 

Four-D nav must be flexible with regard to arrival time. It is envisioned, 
that during a flight the crew would receive periodic updates of the assigned 
metering fix time. When updated, the aircraft computer would determine the 
best speed and altitude for the remainder of the flight to achieve the assigned 
time in a reasonably fuel efficient fashion. It is recommended that subsequent 
R&D develop total 4-D flight, with arrival time flexibility. 

The experience of this test program also highlighted the need for a very 
close collaborative effort with ATC personnel to plan and conduct meaningful 
flight test programs and to develop 4-D techniques that integrate well with 
the ground control function. 


2.5 Future 4-D Flight Management 

The results of the current study can be improved upon with relatively 
simple system improvements. These would include: 

1) Exact aerodynamic/propulsion modeling 

2) Use of RNAV with either radio or radio-inertial position reference 

3) Completely automatic altitude control 

4) Improved vertical damping 

These features will reduce the 2a arrival error dispersion from 19 seconds to 
an estimated 8 seconds. Research by Menga and Erzberger (reference 9) has 
shown that arrival errors of this order are possible. 

The use of direct lift control is another way of maintaining altitude con- 
trol without unduly cycling the engines. It is recommended that subsequent 
research include these fixes in a fuel efficient, flexible 4-D flight management 
approach that combines climb, cruise, and descent from altitude to the runway 
threshold. 


3.0 CONCLUSIONS AND RECOMMENDATIONS 
3.1 Conclusions 

The following are conclusions reached as a result of the activities per- 
formed during this program: 

• The L-1011 demonstrated the capability of executing a 4-D descent and 
achieving 2 a end-of-descent accuracies of 19 seconds with a 2 second 
early mean. 
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• It is estimated that the dispersion can be reduced to 8 seconds with 
a mean of zero seconds utilizing improvements in wind modeling, aero- 
engine modeling, automation of the altitude feedback loop and 
improved procedures. 

• It is also estimated that altitude error at the end of a 4-D descent 
could be less than 50 feet, airspeed error less than 5 kts, and 
position error less than 0.5 n.mi. 

• The 4-D system arrival times should be made flexible so that the end- 
of-descent arrival time can be made to match the ATC assigned metering 
fix time. 

• The 4-D descent appears to decrease air traffic controller workload. 

• Wind modeling should be improved so that arrival time flexibility is 
not severely restricted. 

• A segmented wind model based on actual wind measurements appears to 
offer an improvement over the linear model for stable weather 
conditions . 

• Unstable weather conditions sharply decrease the persistence of wind 
velocity profiles. 


3.2 Recommendations 

The following items are recommended for additional research and develop- 
ment to further the state-of-the-art of 4-D descents: 

• A flexible 4-D descent control law should be designed and flight 
tested to improve ATC integration. 

• The 4-D capability should be extended back through departure, and 
arrival time flexibility worked back as early in the flight as 
possible . 

• Future L-1011 testing should include fully automatic altitude control 
that makes use of both altitude error and vertical velocity feedback. 

• Wind modeling research should be directed toward providing improved 
segmented models. 

• Further analysis should be performed on wind persistence data relating 
to unstable weather conditions. 

• An ongoing dialogue should be maintained with ATC, the airlines, and 
the system developers such as Lockheed and NASA, in order to ensure that 
the progress made satisfies real needs and fits into the ATC system. 
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APPENDIX A 

TIME-IN-DESCENT ALGORITHMS 


This appendix presents the derivation of algorithms describing the 
computation of incremental times required for each of the six segments of 
the L-1011 Flight Management System’s 3-D descent profile. Total time-in 
descent is then simply the sum of each of the six segments. 
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SYMBOL DEFINITIONS 


H 

Aircraft altitude 

(ft) 

H 

c 

Cruise altitude 

(ft) 

h e*d 

E*D altitude 

(ft) 

IAS 

Indicated-Air- Speed 

(kts) 

K 

3600 

(sec/n.mi. /kt) 

K 

c 

Independent Time-Profile "Bending" factor 



(Set equal to 1.0) 


k t 

Temperature Correction Factor 



= (1 + .0043 A Temp) 


R 

Range 

(n.mi. ) 

T 

Time 

(s) 

V 

T 

True-Air-Speed 

(kts) 

v w 

Wind Velocity 

(kts) 

v Wc 

Wind Velocity at Cruise Altitude 

(kts) 

Y 

Flight path angle 

(deg) 

e 

Desired track 

(deg) 

a 

Wind bearing 

(deg) 
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@ E*D SEGMENT 



R, 


R, 


V ar dec 


V h e*d 


V T 

E*D 
IAS 


E*D 


3 

IAS, 

V„ 


l 2 

IAS, 

V 


"1 

IAS. 

V„ 


W, 


W, 


W. 


Level Decel to E*D IAS 

If E*D IAS is equal to descent IAS, set = 0 and omit this step; 
otherwise, compute as follows: 

k(r l ) 


where : 
V„ 


'E*D 


1 ( 1/2 (\* d + V 2 \)j 

= Kj (k 0 +Ki h 1 + k 2 hJ) + (vWV?) IAS 


“ K t (k 0 + K 1 H 1 +K 2 H i) + ( K 3 +K 4 H 1 +K 5 H i) IAS 


E * D j / See Table A-l 
for value of 
constants 
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K 

“r 

K 

c 

IAS. 


3600 sec/n.mi./kt 

1 + .0043 A TEMP 

1.0 


ias 3 - 


(R 3 - Rj 


( IAS 3 


ias e*d) 


IAS 3 = Descent IAS 


W, 


285, 300, 320, 350 or 365, H £ * D > 10000 FT 


250 


, H e * d < 10000 FT 


W 


H 


H Cos ( 0 - or) 


Autopilot Capture 


AT, 


AT, 


where : 
V T 

K 

K m 


k(r 2 - v 


1/2 


'(\*\) “■ v * vV 


K 


k(r 3 - r 2 ) 


1/2 

( V T + V T ) 

1 C os V + V 

& 

> 

+ 


|_\ 1 3 1 

1 



/ 

( 


‘ *T |( K 0 +K 1 H 3 +K 2 H 3) + ( K 3 +K 4 H 3 +K 5 H 3) IAS : 


K 


See Table A- I 


3600 sec/n.mi./kt 
1 + .0043 A TEMP 
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1.0 


K 

c 

IAS 3 


cos v 


= Descent IAS 

10000 FT 
10000 FT 



285, 300, 320, 350 or 365, H „ > 


250 


’ h e*d < 


V r7 - — H, Cos (0 - O') 

H 1 
c 




— — Cos ( 0 - a ) 

H 

c 


(T) CONSTANT 250 KTS IAS SEGMENT 

The equations derived in Segment 4 (Constant IAS) are also applicable to 
this segment with the exception that only the 250 kts IAS coefficients of 
Table A-2 are used. 
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I 



(3) 10,000 FT LEVEL DECEL TO 250 KTS SEGMENT 



AT 


1 


KAR, 


1/2 |V T +V r 


+ 2 \) 


K 

c 
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K AR, 


AT, 


1/2 r( V T 3 + V T 5 ) cos ^ + \ + \' 


■ K 


where: 


\ " *T (K o + 250 V 


T 2 ‘ (K o + 270 V 

T 3 ■ *r (K o + IAS 3 IC 1 ) 


See Table A-4 


= K m 
t 5 T 


|(K o +K 1 H 5 +K 2 H 2 )+(K 3 +K 4 H 5 +K 5 H 2 )IAS 5 1 (See Table A-l) 


IAS 3 = IA£ 5 " (v^J <IAS 5- 270) 


lAS^. = Descent IAS 


= 285, 300, 320, 350 or 365 


j = V = V = V 
W W W W 

1 2 3 4 


(10,000) Cos (0-aO 


V. 


W 


M, 

D 


H 5 Cos (6-«) 


K = 3600 sec/n.mi./kt 


K.^ = 1 + .0043 A Temp 


K =1.0 
c 
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Cos Y = 


R 5 R 3 



+ ( W 


2 


© CONSTANT IAS SEGMENT 


H 

V 


i+1' R i+1' T i+1 

Vi'Xi 


AHj 



K A R. 
x 


AT. = 
x 


1/2 


( V T. + V T. ^ ,) C ° sY i + V + V W. ^ 

\ X X + 1/ 1 X + 1 


where : 
V„ 


T 

i+1 


K 

K 


= K n 


K t (k q + K lH . + K 2 H.) 

( K 0 + K l H i + l + K 2 H i + l) 


3600 sec/n.mi./kt 

1.0 


See Table A-2 


K„ 


= 1 + .0043 A TEMP 




© CONSTANT MACH SEGMENT 
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V w. = H 
l c 

V W 

1 = —± 
W. , . H 
1+1 C 


K = 3600 

*r “ 1 + 


H Cos (0- tv) 

H. , , Cos (8- a) 
l+l 

sec/n.mi. / kt 

.0043 A Temp 


K = 1.0 
c 


Cos V. = 
i 


* R i 



© CRUISE PUSH-OVER SEGMENT 


Above Mach-IAS Cross-Over Altitude 


AT = K 'K +K Mach ) See Table A-5 
c o 1 CRZ | 

Belov Mach-IAS Cross-Ov er Altitude 
AT = K c ''‘ K 0 + ^l IA ^CRZ^ | ^ ee A- 5 

K = 1.0 


c 
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TABLE A-1. V TAS IN E*D ALTITUDE (kts) 


K o 

K 1 

k 2 

K 3 

K 4 

K 5 

3.12522 

■B 

0.26710 x 10' 7 

0.99230 

0.12725 x 1 0' 4 

0.7085 x 1 0' 10 


TABLE A-2. V TAS IN CONSTANT IAS SEGMENT (kts) 


IAS 

K o 

K 1 

k 2 

250 

251.25 

3.318 x 10' 3 

4.525 x 10' 8 

285 

290.00 

3.375 x 10' 3 

6.045 x 10' 8 

300 

301.70 

3.910 x 10‘ 3 

5.300 x IQ' 8 

320 

321.12- 

4.173 x IQ' 3 

5.238 x 10‘ 8 

350 

349.65 

4.715 x 10’ 3 

4.815 x 10' 8 

365 

364.70 

5.011 x 10' 3 

4.710 x IQ' 8 


TABLE A-3. V TAS IN CONSTANT MACH SEGMENT (kts) 


Below 36.090 ft 


Mach 

K o 

K 1 

k 2 

0.800 

525.5 

-1.496 x 10- 3 

-1.004 x 10' 8 

0.820 

542.4 

-1.855 x 10' 3 

-3.912 x IQ' 9 

0.850 

562.5 

-1.923 x 10' 3 

-4.050 x 1 0‘ 9 

1865 

577.2 

-2.306 x 10* 3 

1.908 x 10' 9 


Above 36,090 ft 


Mach 

Vj (kts) 

0.800 

459.0 

0.820 

470.3 

0.850 

485.5 

0.865 

496.0 
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TABLE A-4. V JA$ IN 10 000 ft DECEL TO 250 kts SEGMENT (kts) 


K o 

K 1 

7.81382 

1.12664 


TABLE A-5. AT IN CRUISE PUSH-OVER SEGMENT (seconds) 


Above Mach-IAS Cross-Over Altitude 


Descent Mach 

K o 

K 1 

0.800 

69.85 

-55.17 

0.820 

127.20 

-120.00 

0.850 

165.00 

-164.90 

0.865 

165.00 

-164.90 


Below Mach-IAS Cross-Over Altitude 


Descent IAS (kts) 

K o 

K 1 

285 

156.70 

-0.4292 

300 

156.70 

-0.4292 

320 

130.90 

-0.3258 

350 

114.00 

-0.2400 

365 

95.54 

-0.1559 
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APPENDIX B 


OPEN LOOP DESCENT SIMULATIONS 


This appendix contains computer-generated graphs illustrating the 
effects of descent speed schedule, aircraft weight, temperature and E*D 
altitude on time- in-descent and range for the L-1011. All simulations 
were performed under zero wind conditions. 
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Range from E*D (n. mi) 


- Altitude/range descent profiles, W = 163000 kg 
E*D = 9000 ft/250 kts. 


Figure B-l. 





Altitude (ft x 1(n 








Range from E*D (n. mi.) 


Figure B-5. - Altitude/range descent profiles, W = 163000 kg 
E*D = 11000 ft/250 kts . 
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21 


Figure B-7 


0 120 160 
Range from E*D (n.mi.) 


. - Altitude/range descent profiles, W = 136000 kg, 
E*D = 11000 ft/250 kts . 





Time from E*D (s) 


Figure B-8. - Altitude/time descent profiles, W - 136000 kg, 
E*D = 11000 ft/250 kts . 







0 400 800 1200 1600 2000 

Time from E*D (s) 


Figure B-10. - Altitude/ time descent profiles, W = 163000 kg, 
E*D = 9000 ft/250 kts, M = .82/IAS = 320 kts. 
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APPENDIX C 

CLOSED LOOP DESCENT SIMULATION PLOTS 


This appendix contains computer-generated graphs illustrating the 
results of a closed-loop 4-D descent performed with the wind models of 
Figure 12. The effects upon altitude, descent range and speed are pre- 
sented as well as control surface activity, attitude, and thrust and 
drag requirements. The descent speed schedule used to obtain this data 
was Mach = .82/IAS = 320 kts. 
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B*D (n. mi.) 


rofile, W = 136000 kg, 
) ft/250 kts, Palmdale 










Range from B*D (n. mi.) 


Figure C-3. - Profile descent, W = 136000 kg, E*D - 9000 ft/250 kts, 
Palmdale wind + Tripoli 2<r tailwind. 
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Altitude error (ft) 



9 






Range error (n. mi.) 


I 









Spoiler angle (deg) 
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Normal acceleration (g) 
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APPENDIX D 


L-1011 FLIGHT MANAGEMENT SYSTEM 
DATA BANK INFORMATION 
FOR THE 

DALLAS /FT. WORTH DEMONSTRATION FLIGHTS 


VORTACS 


Name 

I.D. 

Freq. (MHz) 

Latitude 

Longitude 

MVAR 

Hgt (ft) 

Acton 

AQN 

1 10.G 

N32° 

26.1' 

W097° 39.8' 

m 

CO 

o 

848 

Ardmore 

ADM 

116.7 

34 

12.7 

97 10.1 

E9 

930 

Blue Ridge 

BUJ 

114.9 

33 

17.0 

96 21.9 

E8 

610 

Bridgeport 

BPR 

116.5 

33 

14.3 

97 46.0 

E9 

890 

Oailas/Ft. Worth 

DFW 

117.0 

32 

51.9 

97 01.7 

E8 

557 

Gregg County 

GGG 

112.4 

32 

25.1 

94 45.2 

E7 

320 

McAlester 

MLC 

112.0 

34 

51.0 

95 46.9 

E8 

820 

Millsap 

MQP 

117.7 

32 

43.6 

97 59.8 

E9 

890 

Oklahoma City 

OKC 

115.0 

35 

26.5 

97 46.3 

E9 

1360 

Scurry 

SCY 

112.9 

32 

27.9 

96 20.2 

E8 

440 

Shreveport 

SHV 

117.4 

32 

46.3 

93 48.6 

E7 

190 

Texarkana 

TXK 

116.3 

33 

30.8 

94 04.4 

E7 

270 

Tulsa 

TUL 

114.4 

36 

11.8 

95 47.3 

E8 

790 

Waco 

ACT 

115.3 

31 

39.7 

97 

16.1 

E9 

510 

Wichita Fails 

SPS 

112.7 

33 

59.2 

98 35.6 

E10 

1100 

Wink 

INK 

112.1 

31 

52.5 

103 

14.6 

Ell 

2870 
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WAYPOINTS 


I.O. 

Latitude 

Longitude 

HAZEL 

N33° 32.6' 

W095° 18.0' 

ATLAS 

33 25.4 

95 48.7 

(BAT)ON 

33 09.4 

96 36.3 

ALKID 

33 05.4 

96 43.9 


HAMAK 

32 

59.8 

96 

54.4 ( 

BLUE RIOGE FIVE Arrival 

YARBB 

33 

57.7 

96 

07.0 \ 


RADEX 

33 

44.4 

96 

11.8 ) 


HBR 

32 

41.2 

96 

54.4 ' 



Note : HBR is a waypoint in between HAMAK and the approach intercept BATCH for north flow approaches to 
runway 35R. Latitude = same as BATCH; longitude = same as HAMAK 


VELMA 

N34° 

26.6' 

W097° 

30.7' 



MEDIA 

33 

35.2 

97 

39.7 



80IDS 

33 

08.4 

97 

32.3 

f 


PIVIT 

33 

04.5 

97 

23.2 

( 

BOIDS THREE Arrival 

ARINA 

32 

59.6 

97 

11.9 


B22 

33 

32.7 

98 

00.4 



B62 

34 

14.5 

97 

28.0 


1 

AGL 

32 

41.1 

97 

11.9 

/ 

Note: AGL is a waypoint in between ARINA and GEARS for north 

approaches to 35L 

Latitude = 

GEARS; longitude = ARINA 




OONIE 

N31° 

37.1' 

W095° 

43.4' 

V 


OBITS 

32 

03.2 

96 

02.2 



NORMA 

32 

12.3 

95 

05.5 

i 


FINES 

32 

21.3 

95 

46.9 

( 

SCURRY FIVE Arrival 

SEAGO 

32 

37.9 

96 

37.8 


GLADD 

32 

41.1 

96 

43.5 

l 

LANED 

32 

47.3 

96 

54.5 



LGL 

33 

06.2 

96 

54.5 


1 

Note: LGL is a waypoint in between LANED and GARZA for south 

approaches to 1 7L. 

Latitude = 

GARZA; longitude = LANED. 




EDNAS 

N31° 

36.7' 

WQ98 0 

10.1' 

1 

i 

GLENN 

32 

08.0 

97 

50.6 



GRAMY 

32 

24.1 

98 

03.3 

I 

f 

BRYAR 

32 

31.9 

97 

31.8 

( 

, ACTON SIX Arrival 

FLATO 

32 

40.5 

97 

20.1 


CREEK 

32 

46.9 

97 

11.4 

1 

V 

TQA 

32 

14.1 

99 

49.0 



CAR 

33 

06.1 

97 

11.4 


/ 


Note: CAR is a waypoint in between CREEK and ALIGN for south approaches to 17R. 
Latitude = ALIGN; longitude = CREEK. 


WAYPOINTS (cont'd) 


I.D. 

Latitude 

Longitude 

GARZA (intercept) 

N33° 

06.2' 

W097° 

01.8' 

17L (threshold) 

32 

54.6 

97 

01.8 

SL1 

32 

44.6 

97 

01.8 

SL2 

32 

34.6 

97 

01.8 

ALIGN (intercept) 

33 

06.1 

97 

03.0 

1 7 R (threshold) 

32 

54.6 

97 

03.0 

SRI 

32 

44.6 

97 

03.0 

SR2 

32 

34.6 

97 

03.0 

GEARS (intercept) 

32 

41.1 

97 

03.0 

35L (threshold) 

32 

52.7 

97 

03.0 

NL1 

33 

02.7 

97 

03.0 

NL2 

33 

12.7 

97 

03.0 

B(ATC)H (intercept) 

32 

41.2 

97 

01.8 

35R (threshold) 

32 

52.7 

97 

01.8 

NR1 

33 

02.7 

97 

01.8 

NR2 

33 

12.7 

97 

01.8 


Approach to runway 1 7L 


Approach to runway 17R 


Approach to runway 35 L 


Approach to runway 35R 
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WPTS 

EDNAS 

GLENN 

INK 

TQA 

GRAMY 

AQN 

BRYAR 

PLATO 


TUL 

MLC 

YARBB 

RADEX 

TXK 

HAZEL 

ATLAS 

BUJ 

(BAT)ON 

ALKID 

HAMAK 

Tjkc 

VELMA 

B62 

MEDIA 

SPS 

B22 

BPR 

BOIDS 

PIVIT 

ARINA 

SHV 

GGG 

NORMA 

FINES 

DONIE 

OBITS 

SCY 

SEAGO 

GLADD 

LANED 
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ROUTES (cont'd) 


Name I.D. WPTS 


35R, North Flow RNF HAMAK 

HBR 

B(ATC)H 


35R 

NR1 

NR2 


35L, North Flow 

LNF 

ARINA 

AGL 

GEARS 

35L 

NL1 

NL2 

17R, South Flow 

RSF 

CREEK 



CAR 

1 


ALIGN 



17R 



SRI 



SR2 

17L, South Flow 

LSF 

LANED 


LGL 

GARZA 

17L 

SL1 

SL2 
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ROUTE B03 




ARINA 




NL2 A 

I 

I 

NL1 A 

I 



35L 


I 

I 

"A 

GEARS 


A 

I 

A 


NR2 


NR1 


1 

1 


7 

1 

1 

1 


1 

1 

W 


1 


ROUTE BUS 


4 


HAMAK 


A 35R 

I ! I 

A A HBR / 

BATCH _ / 

ROUTE RNF 


DFW NORTH FLOW 


ROUTE RSF 


ROUTE LSF 



CREEK 


A 


ROUTE AQ6 # 


ALIGN 

— A 

I 

17R A 

I 
I 
I 

] 

7 

I 

I 

SRI A 

i 

SR2 A 


GARZA 

A_ 

I 

‘jr 

i 
i 


1 

i 

A SL1 

I 

■ 

SL2 



■ 

A. 


LANED 


ROUTE SC5 


DFW SOUTH FLOW 
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16. Abstract 


This report documents the analysis performed and flight test results obtained by the 
Lockheed-California Company in investigating the feasibility of a 4-D (time-controlled) 
high profile descent. 


Work began in October 1978, to develop a technique to add the time constraint to the 
automatic descent feature of the existing L-1011 aircraft Flight Management System 
(FMS) . Software modifications were incorporated in the FMS computer program and the 
results checked by lab simulation and on a series of eleven test flights culminating 
in a demonstration flight to Dallas/Ft. Worth Regional Airport on August 1, 1979. 

An arrival time dispersion ( 20) of 19 seconds was achieved. Recommended refinements 
should reduce the dispersion to 8 seconds. The 4-D descent technique can be inte- 
grated with the time-based metering method of air traffic control which will be 
introduced in air route traffic control centers in 1980. This combination promises 
substantial reductions in delays at today's busy airports. 
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